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A  NEW  ASYMPTOTIC  THEORY  FOR  THE 
PERIODICALLY  FORCED  LASER 


We  consider  the  laser  rate  equations  for  a  single  mode,  homogeneously 
broadened  and  periodically  modulated  laser.  They  consist  of  two  ordinary 
differencial  equations  for  the  intensity  I  and  the  inversion  of  population 
D.  We  investigate  the  effect  of  periodic  modulations  of  the  cavity  decay 
rate  of  the  form  k(T')  -  «q(1  +  Acos(oT'))  where  A  and  a  are  the  amplitude 


and  the  frequency  of  the  periodic  modulations,  respectively.  In  terms  of 
dimensionless  variables,  the  problem  is  then  formulated  as  [1]: 

-  2I[-1  +  AD  -  Acos(wT) ]  (1) 

^  -  7(1  -  D(1  +  I)]  (2) 

where  T  is  defined  as  T  -  kqT' .  A  is  the  pump  parameter,  7  is  the  loss  rate 


for  the  population  divided  by  the  cavity  decay  rate  and  w  -  o//Cq. 


Our  objective  is  to  analyze  the  periodic  solutions  of  the  laser 
equations  in  the  limit  7  -*  0  (A  >  1) .  This  limit  applies  for  class  B  lasers 

which  includes  ruby,  YAG,  COj  and  semiconductors  lasers  [2].  If  A  >  1  and  A 

-  0,  the  zero  intensity  solution  is  unstable  and  Eqs.  (1)  and  (2)  admit  a 
non-zero  intensity  solution  given  by 


Manuscript  approved  December  24,  1990. 


1 


(3) 


Iq  -  A  -  1  and  Dq  -  1/A. 

A  linear  stability  analysis  of  this  steady  state  solution  shows  that  it 

corresponds  to  a  stable  focus  as  7  -*  0.  In  this  limit,  small  perturbations 

-1  -1/2 
from  the  steady  state  are  decaying  on  an  0(y  )  time  scale  with  0(y  '  ) 

period  oscillations.  This  suggests  to  reformulate  Eqs  (1)  and  (2)  in  terms 

of  the  deviations  from  the  steady  state  solution  (3)  and  to  introduce  a  new 

1/2 

basic  time  t  -  7  T.  These  new  equations  are  formulated  in  [1,3]  and  are 
given  by 

jjj  -  -  y  -  cx  -  e2xy  (4) 

^  -  (1  +  y)[x  -  5cos(fit) ]  (5) 

where  the  new  variables  x,  y,  t,  and  the  parameters  t,  e2>  5,  0  are  defined 
by 

x  -  (D  -  Dq)A<^)1/2,  y  -  (I  -  I0)/I0,  t  -  (27I0)1/2T  (6) 


and 


7l 


e  -  A(rJ-)1/2,  e 2  -  S  -  A(-2-)1/2,  0  -  u>/(27I0) 


1/2 


(7) 


Note  that  «  and  «2  are  0(71/^2)  small  quantities.  The  parameters  S  and  0  are 


two  control  parameters  which  we  assume  0(1). 

Numerical  studies  of  Eqs.  (4)  and  (5)  are  given  in  [3,4],  These 
studies  have  shown  that  the  bifurcation  diagram  of  the  periodically  forced 
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laser  admits  several  branches  of  periodic  solutions  appearing  either  from 
period  doubling  bifurcations  or  from  saddle-node  bifurcation  points. 


THE  CONSERVATIVE  SYSTEM 

If  6  -  «  -  €j  -  0,  Eqs.  (4)  and  (5)  reduces  to  a  conservative  system  of 


equations  given  by 


& 

dt 


-  Y 


£L 

dt 


X(1  +  Y) 


(8) 

(9) 


This  system  is  conservative  and  admits  a  one -parameter  family  of  periodic 
solutions.  A  first  integral  is  given  by 

C  -  Y  -  ln(l  +  Y)  +  hi2  (10) 

where  C  is  the  constant  of  integration.  For  each  C  >  0,  there  exists  a 
periodic  orbit  in  the  phase  plane  (X,Y) .  The  orbit  is  surrounding  the  origin 
and  is  bounded  below  by  the  line  Y  -  -1.  An  example  is  given  in  Figure  la 
for  C  -  4.  The  broken  line  corresponds  to  Y  -  -1.  Figure  lb  shows  Y(t) .  Note 
that  the  oscillations  of  X(t)  are  almost  triangular  and  Y(t)  is  nearly  equal 
to  Y  -  -1  except  during  a  short  interval  of  time.  These  properties  of  the 
periodic  solutions  becomes  more  dramatic  as  C  becomes  larger  and  will  be 
used  in  our  asymptotic  analysis. 


SMALL  MSHBSL  AW  SHALL  AMTUTWE  BBBfilBS 

In  [1]  and  [3],  we  proposed  an  asymptotic  analysis  of  Eqs.  (1)  and  (2) 
which  is  based  on  the  simultaneous  limit  of  small  damping  constants  (e  ■*  0, 
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<2  “  0(e))  and  small  amplitude  forcing  (6  -  0(e)).  In  this  limit,  2im/n- 

periodic  solutions  (n-1,2,3,..)  of  Eq.  (4)  and  (5)  are  constructed  using  a 
regular  perturbation  analysis  in  e .  We  have  found  that  there  exist  for  each 
n  two  distinct  branches  of  periodic  solutions  which  appear  as  S  surpasses  a 
critical  value  S  -  5n(e).  The  branching  of  these  solutions  from  the  limit 

point  located  at  6  -  5^  is  called  primary  saddle  node  bifurcation  in  [3]. 

The  branch  of  saddles  emerging  from  this  point  is  responsible  for  the 
irregular  spiking  observed  in  [2].  The  asymptotic  analysis  based  on  the 
limit  e  -  O(e^)  -  0  is  valid  provided  that  the  values  S  -  $n(«)  are  0(e) 

quantities.  However,  the  exact  numerical  determination  of  the  periodic 
solutions  for  specific  small  values  of  e  and  £2  [3]  indicate  that  the  values 

6  -  is  0(1)  as  n  -  2,  3,..  and  increases  as  n  ■+  <*>.  Moreover,  the  2xn/n- 

periodic  solutions  are  large  amplitude  relaxation  oscillations  similar  to 
the  periodic  solution  of  the  conservative  system  (8)  and  (9)  when  C  is 
large.  We  have  noted  that  these  relaxation  oscillations  already  appear  for 
the  case  n  -  1  when  S  «  0(1) .  These  observations  suggest  to  analyze  the 
periodic  solutions  of  Eqs.  (4)  and  (5)  as  relaxation  oscillations.  To  this 
end,  we  propose  a  new  asymptotic  method  based  on  the  limit  C  -*  <®  keeping  e, 
and  6  fixed. 

THE  PERIODICALLY  MODULATED  CONSERVATIVE  SYSTEM 

In  order  to  clearly  differentiate  the  effects  produced  by  the  periodic 
modulations  and  the  effects  related  to  the  damping  coefficients  e  and  we 
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consider  Che  simplified  problem  with  e  -  -  0.  From  Eqs.  (4)  and  (5),  we 

then  find 


&  -  -  y 

dt  y 


(ID 


^  -  (1  +  y)[x  -  Scos(Ot)] 


(12) 


The  bifurcation  diagram  of  the  2jr/0-periodic  solutions  (period  1  branch)  and 
the  4ir/fl- periodic  solutions  (period  2  branch)  has  been  analyzed  numerically 
for  Q  —  0.9  using  a  continuation  method.  We  have  found  that  the  bifurcation 
diagram  for  e-«2-0isa  good  approximation  of  the  bifurcation  diagram 

when  e  **  0  and  *  0  [3]  if  5  is  not  too  small.  The  period  1  branch  starts 


at  the  origin  of  the  amplitude  vs  S  bifurcation  diagram  and  is  S-shaped.  Its 
left  limit  point  is  located  at  6  -  0.  The  bifurcation  diagram  can  be 
determined  by  a  perturbation  analysis  treating  Eqs.  (11)  and  (12)  as  a 
weakly  perturbed  harmonic  oscillator.  The  method  is  described  in  [1]  and 
correctly  predicts  the  S-shaped  bifurcation  diagram.  The  period  2  branch 
emerges  from  the  period  1  branch  at  a  period  doubling  bifurcation  point  ( S  ~ 
1.26),  is  S-shaped  and  admits  a  right  and  a  left  limit  point  at  6  «  3.61  and 
6-0,  respectively. 

From  the  linearized  theory  for  the  27m/n-periodic  solutions,  it  can  be 
shown  that  the  periodic  solutions  are  either  neutrally  stable  (elliptic 
point)  or  unstable  (saddle  point).  This  suggests  that  Eqs.  (11)  and  (12)  is 
a  conservative  system  of  equations. 
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We  now  propose  to  construct  the  periodic  solutions  using  a  singular 
perturbation  method.  As  C  •*  «,  the  closed  orbit  in  Figure  la  becomes  more 
triangular  and  is  characterized  by  a  slow  evolution  near  y  -  -1  and  a  quick 
large  0(C)  pulse  in  y.  This  suggests  to  construct  the  periodic  solution  in 
two  parts.  The  first  part  (called  outer  solution)  corresponds  to  the  slow 
evolution  of  the  solution.  The  second  part  (called  inner  solution)  describes 
the  rapid  change  of  both  x  and  y. 

Outer  solution.  The  outer  solution  is  characterized  by  y  =  -1.  In  first 
approximation,  we  obtain  the  following  equations 

jg  -  1  and  jg  -  (1  +  y)[x  -  Scos(Ot)].  (13) 

We  solve  these  equations  with  the  initial  conditions 

X(V  “  x0  <  °’  y(V  “  y0  >  *1*  (14) 

After  integration,  we  obtain 

x(t)  -  xQ  +  (t  -  tQ)  (15) 

and 

y(t)  -  -1  +  (y0  +  l)ef(C>  (16) 

where 

f(t)  -  xQ ( t  -  tQ)  +  £(t  -  tQ)2  -  *[sin(0t)  -  sin(0t0)].  (17) 

Since  Xq  <  0,  f(t)  is  negative  near  t  -  tQ  and  y  +  1  is  exponentially  small 
until  f(t)  becomes  positive.  The  critical  time  t  -  t^  satisfies  the 
condition  f(t)  -  0  or  equivalently 

X0(tl  '  t0)  +  2(C1  '  C0)2  '  q! sin(Qci>  ’  sin(0tQ) ]  -  0.  (18) 
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(19) 


At  t  -  t^,  x(t)  and  y(c)  admit  the  following  values 

*<«!,)  “  -  xQ  +  (tx  -  tQ)  and  y(tx)  -  yQ. 

However,  as  soon  as  t  >  t^ ,  y(t)  is  increasing  exponentially  and  the 

approximation  y  »  *1  is  no  more  valid.  Ve  now  assume  that  x  and  y  are 
changing  rapidly  and  propose  a  new  approximation  of  the  solution  valid  for  t 

-  V 

Inner  solution.  The  inner  solution  or  inner  layer  solution  describes  the 
sudden  increase  and  decrease  of  y  and  is  characterized  by  the  fact  that  t 
remains  close  to  t^.  Thus,  this  solution  satisfies  in  first  approximation 

Eqs .  (11)  and  (12)  with  t  -  t^: 


dt  y 


(20) 

(21) 


a*  -  (l  +  y)[x  -  fcosCOe^j. 

We  solve  these  equations  in  the  phase  plane.  A  first  integral  is  given  by 

C  -  y  -  ln(l  +  y)  +  ^x^  -  5xcos(Ot^)  (22) 

where  C  is  the  constant  of  integration.  The  assumption  that  the  inner 
solution  is  characterized  by  large  values  of  both  x  and  y  implies  that  C  is 
large.  We  connect  C  to  the  value  x  •  x^  previously  computed  from  the  outer 

solution  by  using  matching  conditions.  To  this  end,  we  rewrite  (22)  as 

y  +  1  -  expf^^  -  5xcos(0t^)  -  C  +  y]  (23) 

Matching  implies  that  as  y  -*  -1,  x  -*  x^.  Since  C  is  large,  the  condition  for 
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a  bounded  soluCion  requires  Chat  the  expression  in  brackets  is  zero  at  y  -  - 
1  and  x  -  x^.  Thus,  x  -  x^  is  the  positive  root  of  the  following  quadratic 

equation 

2^  -  Sxcos(Ot^)  -  (C  +  1)  «*  0.  (24) 

The  variable  y  is  rapidly  increasing  and  then  decreasing.  From  (21)  and  then 
from  (22),  we  note  that  the  maximum  value  of  y  appears  at  x  -  Scos(flt^)  and 

is  an  0(C)  quantity.  On  the  other  hand,  the  variable  x  is  progressively 
decreasing  from  x  -  x^  >  0  to  x  ■  Xj  <  0  where  y  approaches  y  «  -1.  The 

critical  value  x  -  x^  now  corresponds  to  the  negative  root  of  Eq.  (24)  and 

is  related  to  x^  by  the  following  relation 

x2  -  -  x^  +  25cos(fit1).  (25) 


THE  POINCARE  MAP 

In  summary,  we  have  found  a  map  describing  the  periodic  solution  of  the 
forced  laser  equations  (11)  and  (12).  This  map  connects  the  successive 
minimum  or  maximum  values  of  x.  We  denote  by  x  -  x2n  <  0  and  x  -  x2n+^  >  0, 

the  minimum  and  maximum  values  of  x  corresponding  to  the  times  t  -  t^  and 

tn+^,  respectively  (n-0 , 1 , 2 , . . ) .  From  (18),  (19)  and  (25),  we  have  the 

following  relations  between  x2n>  x2n+l’  Cn  an<*  Cn+1: 


(1)  x2n+l  "  x2n  +  (Cn+l  *  V 


(26) 


<2)  x2n(Vl  '  cn>  *  2(tn+l  '  't/'  nIsin<0Vl)  '  sin«V>  *  0  <27> 


(3)  x2n+2  ‘  '  x2n+l  *  2Scos(nVl) ' 


(28) 
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We  illustrate  our  analysis  by  seeking  particular  solutions  of  Eqs .  (26)- 
(29) .  We  consider  the  cases  of  period  1  solutions  and  the  period  doubling 
bifurcation. 


Period  1  solutions.  Period  1  solutions  satisfy  the  conditions 


x_  .  -  x^  and  t  , .  -  t  - 
2n+2  2n  n+1  n  Q 


(29) 


Inserting  (29)  into  Eqs.  (26) -(28)  leads  to  the  following  relations  between 


x.  ,  x.  ,,  and  t  : 
2n  2n+l  n 


x2n+l  *  *2n  *  <£*>  ' 


x  (^)  +  ^(^)^ 
2n^O  ;  2 ; 


0, 


X2n  "  '  X2n+1  +  2*COS(ntn> 


From  (32)  and  then  from  (31),  we  find  X£n  and  x2n+i: 


(30) 

(31) 

(32) 


x2n  “  ‘  n’  X2n+1  “  fl 


Then  from  (32) ,  we  obtain 


cos(flt  )  -  0 
n 


which  implies 


^"'20 


2nx 
0  ' 


(33) 


(34) 


(35) 


4ir 

x„  .  -  x„  and  t  -  -  t  -  r- . 
2n+4  2n  n+2  n  fl 


The  period  doubling  bifurcation.  Period  2  solutions  satisfy  the  conditions 

(36) 

We  introduce  the  first  condition  in  Eqs.  (26) -(28)  and  formulate  the  problem 
for  X2n>  x2n+i*  x2n+3’  Cn+1  and  Cn+2 '  7116  anatysis  these  equations  can  be 

simplified  if  we  introduce  the  variable  r  defined  as:  r  -  tn+^  "  Cn  '  ^ " 
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r  -  0,  Che  period  2  solution  becomes  a  period  1  solution.  In  terms  of  r,  the 
problem  can  be  reduced  to  two  equations  for  r  and  t^.  Since  the  period 


doubling  bifurcation  occurs  as  r  -*  0,  we  have  analyzed  these  equations  in 
this  limit  and  have  found  the  conditions 


1  +  6fisin(flt  )  -  0  and  cos (fit  )  -  0. 
n  n 


These  conditions  admit  the  solution 


Cn  '  '  Z5'+  ‘  '  l'°- 


(37) 


(38) 


For  Q  -  0.9,  the  period  doubling  bifurcation  is  located  at  S  -  I/O  -  1.1. 
The  exact  numerical  value  is  S  ~  1.26. 


CONCLUSIONS 

In  this  preliminary  report,  we  have  shown  how  to  use  a  singular 
perturbation  approach  to  obtain  a  discrete  map  for  the  driven  conservative 
laser  problem.  Periodic  solutions  of  period  1  and  2  are  located  and  the 
position  of  the  first  period  doubling  bifurcation  agrees  with  the  position 
obtained  numerically  for  the  laser  rate  equations. 

It  is  known  that  the  manifolds  of  the  saddle  periodic  orbit  governs  the 
global  behavior  of  phase  space  such  as  chaos  and  chaotic  bursting.  Since  our 
technique  may  determine  periodic  orbits  which  are  close  to  these  branches  of 
saddle  orbits,  our  laser  map  captures  some  of  the  global  behavior  observed 
in  the  full  system.  The  singular  perturbation  technique  used  here  will  be 
applied  to  the  damped  laser  problem  to  exploit  other  global  features  in  a 
future  paper. 
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